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Abstract 

A mixed II- VI semiconductor Zni_ x Be x Se possesses non-trivial vi- 
bration properties, because its two constituent compounds, ZnSe and 
BeSe, show very different degree of covalency and hence high elastic 
contrast. An anomalous Be-Se vibration line has been observed mostly 
at intermediate Be content in the Raman spectra of thin (Zn,Be)Se 
films. In order to explain microscopic origins and the detailed compo- 
sition of these lines, a first-principles calculation of vibration frequen- 
cies in a mixed crystal has been done, with frozen-phonon technique 
and supercell setup within the density functional theory, by the Siesta 
method, which uses norm-conserving pseudopotentials and strictly lo- 
calized numerical basis functions. The calculations confirmed an earlier 
assumption that the anomalous Be-Se line appears due to the formation 
of continuous chains of a more rigid Be-rich pseudo-continuous phase 
formed within the more soft Zn-rich host region on crossing the Be-Se 
bond percolation threshold (x ~ 0.19). Different local deformation in 
percolated and non-percolated regions affect interatomic elastic inter- 
actions and split corresponding vibration lines. Besides confirming the 
percolation model qualitatively, the calculation provides details about 
vibration patterns in different phonon modes. 
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1. INTRODUCTION 



Chalcogenides of Be, characterized by a strong covalent character of their 
chemical bonding, are outstanding among the II- VI semiconductors which are 
otherwise largely ionic systems. BeSe makes solid solution with ZnSe in the 
whole composition region. Beyond providing the possibi lity to tune the funda- 
mental band gap and dielectric function in broad limits (jWilmers et al. 103), 
the coexistence of strongly ionic and strongly covalent bonds in the same crys- 
tal makes its elastic and vibrational properties interesting. The frequencies of 
optical modes in pure ZnSe and BeSe are far separated. This remains basi- 



cally true in the alloy ([Pages et all l200nh . The remarkable point is that an 



additional Be-Se mode has been detected, mostly in the intermediate concen- 
tration range, 20 t o 80 at.% Be. This effect has been tentatively attributed by 
Pages et all (j200l| ) to the onset of percolation on the Be sites, i.e., creation of 
infinite -(Be-Se)- chains in the predominantly ZnSe crystal. 

In the attempts to understand and describe the lattice dynamics of mixed 
crystals in terms of simple models, an obvious difficulty is that interaction 
parameters (force constants) between two given atoms depend on their local 
environment in a not straightforward way. The elastic couplings are affected 
by interatomic distances which may differ in an alloy; moreover the combi- 
natorial aspect (how many neighbours of one or another type) may play an 
unknown role. Therefore it could be difficult to agree on a simple and yet 
realistic mode l. First-principles ca l culations of l a ttice vibrations, either by lin- 



ear response ( Pavone et a/. , 19931 Qonze et al . 1994 : Baroni et al . 2001 ) or 



frozen-phonon ([W ei an d Cho 



(Parlinski and Kawazoe 



schemes, work well for perfect crystals 



1999) or ordered mixed phases ( Sternik et al . 19991 : 



Baroni et all l200lh but face problems in case of substitutional disorder. In 



this latter case, the treatment of representative supercells remains a virtually 
unique alternative, with the obvious shortcoming that the selection of such su- 
percells is limited, and the calculation technically demanding. In the present 
contribution, we calculate electronic structure, force constants and phonon fre- 
quencies in the (Zn,Be)Se mixed system usi ng the first-principles method, and 
calculation code, Siesta hw^l B ^ calculation; were practically 
done for Zn 3 BeSe 4 , ZnBe 3 Se 4 and Zn26Be 6 Se 32 supercells, the latter simulating 
both percolating and "isolated" Be ions in the lattice. 

The additional vibrational modes, whose appearance in the calculation 
could be traced back to the onset of percolate d Be chains, closely resemble 
those observed earlier in experimental spectra bv lPages et al. (|2002b|) and out- 
lined briefly in Sec. 2. This supports earlier assumptions about the role of per- 
colation in the formation of split-off vibration lines. The results of calculation, 
which were organized as described in Sec. 3, allow to study the composition of 
each vibration mode separately, and to explain the underlying microscopies, in 
terms of variation of interatomic distances. The discussion of results is given 
in Sec. 4. 
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2. EXPERIMENTAL SITUATION 



The vibrations spectra of (Zn,Be)Se alloys have been detected by Raman mea- 
surements in the backscattering geometry from thin films grown by molecu- 
lar beam epitaxy on (001) GaAs. The experimental setu p, obtained spectra 



and their discussion have been in detail described earlier (jPages et all |2 



Pages et al . 2001 . l2002a| ). A typical selection of spectra, in dependence on 
concentrations, is shown in Fig. [TJ For the sake of our subsequent discussion, 
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Figure 1: Raman spectra of Zni-^Be^Se films grown on the GaAs substrate, mea- 
sured in the backscattering geometry along the (110) edge axis. The "regular" and 
"anomalous" (l abeled TO| g ) l ines coexist in the concentration range 20 to 80 
at.% of Be. See lPages efnU {2001) for more detailed discussion of the experiment. 
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we emphasize the region of transverse optical (TO) Be-Se modes, in the range 
400-500 cm -1 . The Zn-Se mode, around 200 cm -1 , is useless beca use of strong 



Fano interference with a disorder-activated acoustical continuum (Pages et al 



2001). Incidentally the large frequency gap between the Be-Se and Zn-Se 



nodes is basically attributed to the large difference in the reduced mass of two 
bonds, in the ratio 1:4. The increase in the Be concentration shifts the BeSe- 
like TO line upwards. From ~20 % of Be on, an additional ( "anomalous" ) peak 
is formed at ~70 cm -1 on the low-frequency side from the main one. Two peaks 
coexist through the concentration range ~20 to ~80% of Be, both displacing 
towards higher frequencies and approximately maintaining their separation. 
The intensity of the anomalous peak grows at the expense of the normal one, 
until finally (on the Be-rich side) the anomalous peak remains the only one, 
and ends up as the narrow TO line in the limit of pure BeSe. 

The discussion in previous works ([Pages et all l200lL l2002bf) emphasized 
that such behaviour cannot be accounted for by segregation of crystal into 
BeSe-rich and ZnSe-rich fractions. The proposed explanation of the anomaly 
was the following. Above the Be-Se bond percolation threshold (~19 at.% Be), 
a pseudo-infinite wall-to-wall chain of Be-Se bonds is formed spontaneously 
throughout the alloy. Based on simple local strain considerations, the bond 
length within the chain was presumed to be larger than the bond length of the 
Be-Se bonds isolated in the ZnSe-like host matrix, with concomitant impact 
upon the phonon frequencies. However, microscopic justification of this picture 
was so far lacking. Our present aim is to construct a representative structure 
model, compare the vibration properties following from it with experimental 
observations, and in case of agreement to analyze microscopical origins of the 
observed anomaly. 



3. CALCULATION METHOD AND SETUP 



First-principles calculations have been performed wi thin the general scope of 
the density functional theo ry (see, e.g. , Kohn . 19991 for a review), using the 
calculation method Siesta ( Soler et al. . 2002f) with its compact basis of strictly 
confi ned atom-centered functions ((Sanchez- Portal et al . 19961: Junquera et al . 
2001). Double-C basis functions with polarization orbitals were used for all 



valence states - Be 2s and 2p, along with 3d, 4s and 4p for Zn and Se. Note 
that we attributed the semicore Se3d states to the val ence band. The norm- 
conser ving pseudopotentials have been generated after Troullier and Martins! 
( 19911 ). The results discussed below have been obtained using the local density 
approximation for the exchange-correlation. 

The calculations have been done for three different cubic supercells, shown 
in Fig. El The smaller ones included four formula units in the zincblende 
structure. BeZn3Se4, with the nominal Be concentration (25 at.%) already 
beyond the percolation threshold, contained nevertheless only isolated (and 
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Figure 2: BeZn 3 Se4, Be3ZnSe4 and BeeZn 2 6Se32 supercells used in the calculation. 
More covalent Be-Se bonds are shown by thick lines. The large supercell contains a 
continuous (-Be-Se-) chain and two isolated Be atoms (indicated by arrows), who 
have only Zn atoms as next-nearest neighbours. 

artificially ordered) Be impurities, whose second-nearest neighbours are Zn 
atoms only. The "symmetric" case of 75 at.% Be (BesZnSe^ contains, on the 
contrary, isolated Zn atoms and fully connected Be-Se network. Moreover, 
a larger supercell Be 6 Zn 2 6Se 3 2 simulated the onset of percolation on the Be 
sublattice, i.e., the formation of continuous (-Be-Se-) chains, along with the 
presence of isolated Be substitutions. The Be concentration in this case is 
just above the percolation threshold on the zincblende lattice. In all cases, 
the calculated forces and strain were used for relaxing the lattice vectors and 
internal coordinates by the conjugated gradient algorithm, without any sym- 
metry constraints imposed. This allowed us to discuss the variation of the 
lattice parameters and interatomic distances in their dependence on concen- 
tration and local environment. The vibration frequencies have been calculated 
in the frozen phonon approach, sequentially introducing finite displacements 
of atoms in the supercell by 0.03 Bohr along three Cartesian directions from 
their equilibrium positions, and analyzing the forces induced on all atoms. In 
a single cell, this technique gives access to the zone-center TO phonon only, 
but in a larger supercell the analysis of phonon eigenvectors allows to make 
some conclusions about the phonon dispersion. 

As benchmarks for testing the quality of the pseudopotentials used and 
of the basis sets, we calculated the equilibrium lattice parameter and the T- 
TO frequency in pure ZnSe and BeSe. The calculation results are shown 
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Table 1: Equilibrium lattice constant a, bulk modulus B, and frequency of the T- 
TO phonon uj for pure ZnSe and BeSe from the present calculations (WIEN97 and 
Siesta), previous calculations, and experiment. 



Method 




ZnSe 






BeSe 




a (A) 


B (Kbar) 


uj (cm x ) 


a (A) 


B (Kbar) 


uj (cm 1 ) 


WIEN97 


5.571 


727 


162 


5.087 


831 


456 


Siesta 


5.590 


778 


200 


5.114 


965 


429 


exp. 


5.668 


624-647 


207" 


5.137 b 


920 


501 c 


other calc 




750 d 


216 e 









ISzvbowicz eTal\ (l200r#) : iLuo et all lll995lk IPages et all <l200(t : ICohenl Jl98.^ : 
faennion efal\ l|l97l|) 

in Table I along with the experimental data and the results of another, all- 
electron calculation, which has been performed using an accurat e all-electron 



full-p otential linearized augmented plane wave method WIEN97 fBl aha et al. 



1997). The difference between the Siesta results and the experiment reflects 
the shortcomings of the local density approximation, and moreover technical 
limitations of the pseudopotential generation and the choice of basis functions. 
A similarly more faithful estimation of some results with Siesta, as compared 
to WIEN97, in view of numerical superiority of the latter method, is probably 
due to some error cancellation. Anyway, Table I gives an idea of the level 
of accuracy of our present calculation approach. It is obvious already from 
comparison the bulk moduli that ZnSe is much softer than BeSe. The same 
trend holds for the force constants, and mainly accounts for the difference in 
the T-TO phonon frequencies. The difference in Be and Zn masses further 
enlarges this difference. 



4. RESULTS AND DISCUSSION 

Fig. El shows the calculated equilibrium lattice constant (derived from the vol- 
ume of relaxed supercells) and nearest-neighbour distances in pure materials 
and three mixed- composition supercells. One can see that the (Zn,Be)Se lat- 
tice parameter quite accurately follows the Vegard's law. At the same time, 
the Be-Se and Zn-Se bond lengths, which differ by about 9% in the parent 
compounds, tend to remain diverse, and almost independent on concentration. 
It means that the accommodation of short Be-Se and long Zn-Se bonds in the 
lattice of a mixed crystal may lead to local strains and frustrations. In fact, in 
the large supercell one finds a distribution of bond lengths of one and another 
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relaxed structure of (Be,Zn)Se supercells 




Figure 3: Mean lattice constant (top panel) and nearest-neighbor bond lengths 
(bottom panel), according to calculations on two pure end systems and three mixed- 
composition supercells. In the BeeZn26Se32 supercell, the distribution of Be-Se and 
Zn-Se bond lengths is shown. 

type, shown as continuous curves in the bottom panel of Fig. EJ 1 A further 
analysis of these different bond lengths shows that isolated (non-chain) Be 
atoms can easily induce a local and "isotropic" lattice contraction, binding all 
their four Se neighbours at a distance of 2.23 A, i.e. at nearly the nominal 

1 Obviously the number of nearest-neighbour bonds in the 64- atom supercell is finite 
(=128); an artificial smearing of the bond distribution function is introduced for better 
visibility. 
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bond length of pure BeSe. About the same, or even shorter (2.21 to 2.23 A), 
bond length gets stabilized between Be atoms in the chains and their off-chain 
Se neighbours. One finds finally that the Be-Se bonds along the chains are 
the most difficult to shorten, and exhibit a variation of lengths 2.24 - 2.28 A. 
This can be understood because the chain is infinite, and its period is essen- 
tially fixed by the mean lattice constant of a crystal with a predominant Zn 
concentration. The only possibility to shorten the chain links is at the price 
of increasing the Se-Be-Se angles, making a zigzag chain more straight. This 
increases average distances between in-chain Se atoms and their Zn neighbours 
beyond the comfortable ZnSe bond length. 

The discussed differences in interatomic distances, and their stretching be- 
yond, or below, the natural bond lengths has an effect on the calculated force 
constants. A diagonalisation of the dynamical matrix yields frequencies and 
eigenvectors of all 192 zone-center vibration modes of the supercell, includ- 
ing three acoustic modes (the closeness of whose frequencies to zero is a good 
test for the numerical accuracy of the calculated forces). Fig. 0] shows phonon 
densities of states, constructed as 



where Ui are (discrete) vibration frequencies in the supercell (leaving aside the 
three zero-frequency acoustic modes), A a are elements of the eigenvector of 
the corresponding mode, and the summation is done over {a}, a selection of 
atoms in a certain structure category. A broadening parameter A=10 cm -1 
was introduced for better visibility. 

The top panel of Fig. 0] shows the contributions of Zn and Be atoms, which 
are well separated in the frequency domain. The lower p art (up to ~300 cm~* ) 
resembles the phonon density of states in ZnSe (see, e.g., Hennion et al . 19Tlh . 



\0Ji 



The high-frequency part represents a Be-Se vibration; its fine structure reveals 
the differences in the force constants which couple different Be-Se pairs. The 
vibration modes involving the isolated Be are densely situated around 490 
cm -1 . On the contrary, the vibration bands involving Be atoms in the chains 
split into a softer and a harder components. A comparison with the phonon 
DOS projected on Se atoms (Fig. EJ bottom panel) shows that the softer one 
affects almost exclusively the Se atoms in the chains (i.e those with two Be 
and two Zn neighbours), whereas the harder one involves also the peripheric 
Se atoms, which build just one bond to Be in the chains. This separation 
immediately follows from the distinction between extended Be-Se bonds along 
the chain, with correspondingly weakened force constants, and "natural" Be- 
Se bond lengths to off-the-chain Se atoms, or around isolated Be impurities. 
This observation clearly indicates that the splitting off of a softer "additional" 
line from the main TO Be-Se peak in the measured spectra (Fig. [TJ is due to 
the formation of continuous (-Be-Se-) chains, where the Be-Se bond length is 
essentially fixed by the predominantly ZnSe bulk, therefore extended over the 
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Be 6 Zn 26 Se 32 supercell: phonon DOS 




Figure 4: Phonon densities of states for different categories of atoms represented 
in the BeeZn26Se32 supercell. Top panel: Be and Zn contributions, bottom panel: 
Se contributions. Three acoustic modes at co=0 are removed, and the broadening of 
discrete vibration lines is introduced. 

"natural" Be-Se distance, and results in weakened force constant. The drift of 
the "normal" and "additional" lines with concentration, seen in Fig. is not 
directly addressed by our simulation on a single "percolated" supercell, but it 
can be easily understood from the general trend of the average lattice constant. 
A gradual shortening of Be-Se distances with the increase in Be concentration 
hardens all vibration modes. The concept of percolation on the Be sublattice 
also makes clear why the "additional" and not the "normal" mode transforms 
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into the single TO mode in the limit of high Be concentration (see Fig. [Q: the 
"additional" mode was that involving chained (-Be-Se-) structures, which at 
high Be concentration become omnipresent. 

In order to get a better idea of vibration patterns within "normal" and 
"additional" lines, we shown in Fig. El the snapshots of four vibration modes. 
Whereas the most of vibration modes involve many atoms in the supercell in 
a complicated way, these four allow more or less easy understanding of their 
underlying displacement. Such simplified schemas are shown below each mode, 
neglecting the bending of the zigzag chain. The most general observation is 
that the movement of Be atom in the common plane with two Se neighbours 
occurs at lower frequency than out-of-plane vibrations. This is however a quite 
simplifying statement; one should consider that the off-the-chain Se atoms also 
participate in these movements, up to different extent in harder ("normal") 
and sifter ( "additional" ) modes, as has been mentioned above. 

We leave for a subsequent more detailed discussion elsewhere the analysis 
of vibrations in smaller supercells, which we used here only as a reference to 
the variation of lattice parameters. Moreover, the analysis of vibration modes 
with respect to their wave vector is important and will be reported separately, 
as only q=0 vibrations contribute to the Raman spectra. 




Figure 5: Snapshots of four selected vibration modes involving Be atoms in the 
chains. The displacements of Se atoms are negligible in the chosen scale. For each 
mode, a simplified schema of its characteristic displacement is shown below. 
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5. SUMMARY 



By performing a straightforward calculation of vibration frequencies in a large 
supercell simulating (Zn,Be)Se mixed system, we were able to explain the onset 
of "anomalous" vibration mode in the frequency range close to the T-TO line of 
BeSe. This mode appears due to a softening of corresponding force constants 
in expanded Be-Se bonds, that occurs as continuous Be-Se chains start to 
traverse the Be^Zn^Se crystal on reaching the Be percolation limit of x—19 
at.%. We found remarkable variations of the Be-Se bond length in a (Zn,Be)Se 
mixed crystal with the Zn predominance. The Be-Se distances shorten most 
efficiently between Be ions in the percolating chains and their neighbouring Se 
ions not belonging to the chains. On the contrary, the Be-Se distances within 
the chains are essentially fixed by the ZnSe bulk, i.e. they are over-stretched, 
as compared to relaxed BeSe crystal. The Be-Se bond length for a Be ion being 
an isolated substitutional impurity in the ZnSe matrix is intermediate between 
two previously mentioned cases. As a consequence, one finds distinct Be-Se 
force constants for these three cases, and hence the splittings of corresponding 
vibration modes. 
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